Spectral scrambling in Coulomb-blockade quantum dots 

Y. Alhassid 1 and Yuval Gefen 1,2 
1 Center for Theoretical Physics, Sloane Physics Laboratory, Yale University, New Haven, 

Connecticut 06520, USA 
Department of Condensed Matter Physics, The Weizmann Institute of Science, Rehovot 76100, 

Israel 
(August 20, 2001) 

Abstract 

We study the fluctuations of an energy level as a function of the number 
of electrons m added to a Coulomb-blockade quantum dot. A microscopic 
calculation in the limit of Koopmans' theorem predicts that the standard 
deviation of these fluctuations behaves as y/m in the absence of surface charge 
but is linear in m when the effect of a surface charge in a finite geometry 
is included. The microscopic results are compared to a parametric random- 
matrix approach. We estimate the number of electrons it takes to scramble the 
spectrum completely in terms of the interaction strength, the dimensionless 

Thouless conductance, and the symmetry class. 
PACS numbers: 73.23.Hk, 05.45+b, 73.21.La, 73.23.-b 
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The simplest model for describing a quantum dot in the Coulomb-blockade regime is 
the constant interaction (CI) model in which the electrons occupy single-particle levels 
and the Coulomb interaction is taken as an average electrostatic energy that depends only 
on the number of electrons. When the single-particle dynamics in the dot are chaotic (or 
diffusive), one can apply random matrix theory (RMT) to describe the statistical properties 
of the single-particle wave functions and energies within an energy window whose width 
is the Thouless energy. RMT was successful in describing the conductance peak height 
distributions and their sensitivity to time- reversal symmetry However, other measured 
statistics, most notably the peak spacing distribution indicated that it is necessary to 
include interaction effects beyond the simple CI model ||. The best way to take into 
account interactions while retaining a single-particle framework is the Hartree-Fock (HF) 
approximation, which has been used to explain some of the observed features of the peak 
spacing statistics ||. 

The HF single-particle wave functions and energies are calculated self-consistently and 
therefore can change as electrons are added to the dot. In the statistical regime (i.e., in a 
chaotic or diffusive dot), this phenomenon is called scrambling. Scrambling was observed 
in the decay of correlations between the m-th excited state in the dot and the ground-state 
of a dot with an additional m electrons 0, and indirectly through the saturation of the 
peak-to-peak correlations as a function of temperature HH. A phenomenological way to 
describe scrambling is to consider a discrete set of Hamiltonians (corresponding to the differ- 
ent number of electrons in the dot) that are random but have the correct symmetries. Such 
a set is known as a discrete Gaussian process (GP), and can be embedded in a continuous 
GP, i.e., random matrices that depend on a continuous parameter [|K|. This parametric ap- 



proach leads to a nearly Gaussian peak spacing distribution [|TTJ , and explains the saturation 
of the number of correlated peaks versus temperature ||. While the parametric approach 
is appealing in its simplicity, it is not clear how well it describes features obtained in a 
microscopic approach. For example, in the parametric approach, the fluctuation standard 
deviation (FSD) of a typical level upon the addition of m electrons into the dot is propor- 



tional to m. One can also evaluate the FSD of a level in the microscopic approach through 
the fluctuations of the diagonal interaction matrix elements (in the limit where the single- 
particle wave functions do not change with the addition of electrons, hence the off-diagonal 
matrix elements can be neglected). When only the bulk screened Coulomb interaction is ac- 
counted, we find that the associated FSD behaves like \/m, unlike the standard parametric 
approach. However, in the finite geometry of a quantum dot, the level fluctuation may be 
dominated by the variation of the mean-field potential arising from surface charge. We shall 
show that in this case the FSD is linear in m (for 1 < m « j, where g is the dimension- 
less Thouless conductance), in overall agreement with parametric RMT. The main results 
of the microscopic approach are summarized in Eqs. (H), ([12]) and ( |13|) in the absence of 



surface charge and Eqs. fllTf ), (|18|) and ([19|) in the presence of surface charge. Extrapolating 
our results to larger values of m, we also estimate the dependence of the number of added 
electrons m c required for complete scrambling on g and the interaction strength both in the 
absence and presence of surface charge (Eqs. (fH| ) and (pT|), respectively). 

Scrambling implies variation of both eigenstates and energy levels. However, in this work 
we approach the problem from the limit where the wave functions do not change (Koopmans' 
theorem |JT2|), and we only estimate the FSD of energy levels. In the parametric approach 
this limit is valid when the change in the parameter upon the addition of m electrons is small 
compared to the mean parametric distance between avoided crossings. Complete scrambling 
occurs when the energy fluctuations become comparable to A (the mean-level spacing in the 
middle of the band). 

We first discuss the parametric approach. The variation of the single-particle energies 
and eigenfunctions (e.g., in a mean-field approximation) with the addition of electrons into 
the dot is described by a parametric variation of the Hamiltonian We denote by H(xjy) 
the effective single-particle Hamiltonian of the dot with M electrons. We assume that the 
dot is either diffusive or ballistic with chaotic dynamics, and that the statistical properties 
of H(x_v) are not modified by the interactions. Restricting ourselves to the universal regime, 



i.e., to ~ g levels in the vicinity of the Fermi energy [14], we assume that H(xj^) belongs to 



one of the Gaussian ensembles of random matrices whose symmetry class (3 is independent 
of x^. The sequence of Hamiltonians H(xj^) forms a discrete GP that can be embedded 
in a continuous GP H(x) A simple GP is given by [p~Q| , p~9|] H(x) = cosxifi + sina;iJ2, 



where Hi and H2 are N x N uncorrelated random matrices chosen from the Gaussian 



ensemble of symmetry class (3: P{H) oc e~2^ TriT . We choose a = (2/N) 1 ^ 2 ( so that 
the average level density (a semicircle) has a constant band width of 2a\j2N = A( and 
the mean- level spacing in the middle of the spectrum is A = T\aj\/2N = 7r(/N. The 
average distance between avoided crossings is given by the inverse of the rms level velocity 



5x c = A[((9e a /(9x) 2 ] -1 / 2 = tt(P/2N) 1 ^ 2 , where e a is an energy level. This distance is larger 
for the GUE by a factor of \[2 compared with the GOE. Two-point parametric correlators 
become universal when the parameter x is measured in units of 8x c ; 1.6. ; cLS cL function of a 
scaled parameter x = x/Sx c . 

The energy levels e a scramble as the parameter x changes: 5e a = e a (x + Sx) — e a (x). The 
variance of 5e a is estimated in the limit Sx 1 to be a 2 (6e a ) = A 2 (5x) 2 using first order 
perturbation theory in 5x (i.e., ignoring the change of the single-particle wave function as 
x — > x + Sx). 

In the parametric approach it is assumed that x changes by Sxi upon the addition of one 
electron into the dot (independently of the number of electrons J\f). Thus for the addition of 
m electrons Sx m = mSxi, and as long as 5x m 1, we can still use first order perturbation 
theory to find 

criSe^) = AmSxi = ma (6 eg*) , (1) 

where 5e^ = e^ v ' +m ) — (e^ is the energy of level a in a dot with J\f electrons). 

To relate the parametric approach to a microscopic mean-field approach, we describe the 
effect that adding electrons has on a particular single-particle HF level. The Hamiltonian of 
the dot is 

H = 5 h i? a l a 3 + v fjki a Wj^k , (2) 

ij ijkl 



where h^ ' is the one-body part (accounting for the disorder or the chaotic confining po- 
tential) and v^. kl = (ij\v\kl) — (ij\v\lk) are the antisymmetrized matrix elements of the 
two-body interaction. The HF Hamiltonian in the self-consistent basis for M electrons is 
given by hffl = + YjgLi v aS^s (where the sum is over the lowest M occupied levels), and 
the eigenvalues of hffl are the HF single-particle energies e^f\ The diagonal part of the 
Hamiltonian (0) in the HF basis is 

^diagonal = + \Y, V i^ n ^l > ( 3 ) 

a a, 7 

where n a is the number operator for single-particle state a and ud, = . 

In Koopmans' limit, the single-particle wave functions do not change when an electron is 
added to the dot. Consequently, the change in the single-particle energy e a when an electron 
is added to the M + 1-st level is given by a diagonal interaction matrix element 

s (1) = IN+l) _ (AT) ^ A /.N 

Assuming the HF wave functions satisfy similar statistics as in the non-interacting case, 
we can estimate the variance of the matrix element in (HI). The two-body interaction used in 



the HF approximation is often an effective interaction, e.g., an RPA screened interaction |16 



Alternatively, one can employ a short-range dressed interaction [|T7| v(r—r') = XAV5(r—r'), 
where A is the interaction strength and V is the system's volume. A diagonal interaction 
matrix element is then given by v ai = XAV J dril)^(r)i/;*(r)i/; a (r)i/; 7 (r), and its variance 



a 2 (v aj )= X 2 A 2 V 2 J dn J dr 2 [^(rO^^rO^^n)^^)^^)^^)^^^)^^); 

[ra(riW 7 (ri)Mri)^(ri)) (^2)^2)^2)^2))] • (5) 



Only the connected part of the ensemble average over the product of eight wave functions 
contributes to the r.h.s. of (|5]). In the following we restrict our calculations to the GUE 
symmetry. All possible pairwise and quadruplet-wise contractions should be taken into 
account, and they can be regrouped to give 
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a 2 (v ay )= \ 2 A 2 V 2 J dv x j dr 2 { [(C(n)Va(r-i)^a(r 2 )C(^)} {^{rx)^{r x )^{r 2 )^{r^) - 
+ [(C(n)^(r 1 )^ 7 (r 2 )^(r 2 )) (^(r0^ 7 (ri)^(r 2 )C(r 2 )) - ^ 
+ (C(r-i)^(r- 2 )^ 7 (r 1 )^(r 2 )) (^1)^2)^1)^2))} , 

where the l/V 4 factors are subtracted to avoid double counting of pairwise contractions 
(here (. . .) denotes an ensemble average and not just the connected part). 

Wave- function correlations appearing in Eq. @ were calculated in Ref . K2^j2l| • Including 
contributions of order 1/g 2 and assuming that the distance between levels a and 7 is less 
than the Thouless energy, we have 

, A 2 A 2 



o- 2 (v ai )= — 



j dr x J dr 2 {[l + k d (r)(l + xl(r 1 ,r 2 )) + U(r 1 ,r 2 )+b(r 1 ,r 2 )] 2 -l 



+ [1 + k d (r)U(r u r 2 ) + c(r ls r 2 )] 2 - 1 + [k d (r) + n(r ls r 2 )] 2 } . (7) 

The function kd(r) (r = r\ — r 2 ) is given by k d {r) = (7rz/) 2 (ImG i? (r)) 2 (where G R 
is the retarded Green function and v the electron density of states per unit vol- 
ume) and describes the short-range correlations. In 2D k d {r) = exp(— r/£) Jq (k F r) 
(where I is the mean free path and Jo is a Bessel function), while in 3D k d (r) = 
exp(-r/T)(sin k F r/k F r) 2 . Tl(r 1 ,r 2 ) = (iris)' 1 Eq/o 4>Q{ r i) ( t ) Q{ r 2) / DQ 2 is the long-range 
diffuson propagator where <j)Q is the eigenfunction of the diffusion operator correspond- 
ing to the eigenvalue DQ 2 (D is the diffusion constant). The functions b, c in Eq. (|7|) 
are of order 1/g 2 (c ~ U 2 (r 1 ,r 2 )/2 for I <C r < I |2l] where L is the linear 
size of the dot). The main contributions to (0) arise from the long-range terms II 2 
(/ c/r 1 II(r 1 , r 2 ) = and does not contribute). Using the normalization conditions 
Jdr 1 jdr 2 (C(n)^(n)^(r- 2 )C( r 2)> = 1 and Jd ri Jdr 2 (CW^nW^)^^)) = 
1, it follows that J drx I dr 2 [k d {r) (1 + II(ri, r 2 )) + b(r lt r 2 )] = and 

/ dri J dr 2 [1 + k d (r)H(rx, r 2 ) + c(r x , r 2 )] = 0, respectively. Thus the functions b and c 
do not contribute to order 1/g 2 and 

a 2 (v a ,) = 2K(XA/g) 2 , (8) 
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where g = 2itvDL d ~ 2 is the dimensionless Thouless conductance. For a cube in d dimensions 
K = \ Zm^o l/l n l 4 5 where the summation is over n = (rii, . . . ,n d ) with rii non-negative 
integers [n ^ (0, . . . , 0)]. 

We next consider the addition of m electrons into the dot. Extending Koopmans' limit 
to m electrons, i.e., assuming the single-particle wave functions do not change with the 
addition of m electrons, we can express the change in a single-particle level a as 

m 

^ = e +m) -^£^ +l - (9) 

8=1 

For m < g, the variance of (§) is given by a 2 (5e^) = ma 2 {v ai ) +m(m — l)(v^v^s — v^ v^s) 
(since the variances and covariances are independent of the particular wave functions). In 
analogy to Eq. @, the covariance of v ai and v a s is given by 



1 



\ 2 A 2 V 2 J dr.J dr 2 



(C(r-i)^a(r-i)^a(r-2)C(^))<^(^i)V'7(^i)^(r-2)VK^) / yi 



+ 



(VS(ri)^«(ri)^(r a )^(ra)> (V?(ri)^(riW a (r a )<(r a )) - ^ 



(10) 



and using the known wave-function correlations 
A 2 A 2 



J dn J dr 2 {[1 + k d (r)(l + U(r lt r 2 )) + U(r x , r 2 ) + 6(r 1} r 2 )] 



^ 2 

x [1 + ^(r)n(n, r 2 ) + c(n, r 2 )\ - 1 + [1 + fc d (r)n(ri, r 2 ) + c(n, r 2 )] 2 - l} . (11) 

Terms such as fcjjll, kjU 2 and fc^c m Eq. ( |TT| ) contribute ~ (£/L) 2 X 2 A 2 /g 3 and behave like 
~ X 2 A 2 /g 3 in the ballistic limit £ ~ L. Other terms contribute ~ l/p 4 - Thus to leading 
order in 1/g 2 

VU-^tU~0. (12) 
Using Eqs. (H), (H) and QT2) we find E3| 



<x 2 (^ m) ) = ma 2 (5ei 1} ) = 2m/t(AA/( ? ) 2 . (13) 



7 



We conclude that the statistical fluctuations accumulated by sequentially adding electrons 
into the dot are added independently, namely, the FSD of an energy level behaves like y/m, 
contrary to the parametric approach where it is linear in m (see Eq. ([[])). 

A complete scrambling of the spectrum corresponds to a number of electrons m c for 
which a(5e^) ~ A. Using (|TJ) we find 

™ c ~(£f- (14) 
We note that since we have neglected wave-function scrambling due to off-diagonal inter- 



action matrix elements, our estimate in Eq. ([14]) should only be regarded as an upper 
bound. 

Next we consider the effect due to excess negative charge on the boundaries of the dot. 
The effective interaction due to such charge is [H| v(r,r') = V K (r) + V K (r') where V K (r) 
describes the variation of the mean-field potential upon the addition of one electron. For 
a disk of radius R in 2D we have V K (r) = — (e 2 /2kR)(R 2 — r 2 ) -1 / 2 where n is the inverse 
screening length. The related interaction matrix element is now v ai = J dr\ifj a (r)\ 2 V K (r) + 
/ dr\ip^(r)\ 2 V K (r) . Using the wavefunction correlations to order 1/g 2 , its variance is 

° 2 M = TT2 f dr i / dr 2 [(v 2 (\Mri)\ 2 \Mr2)\ 2 ) - l) + (V 2 ^^)^^)?) ~ 1 



xKNKN = ^ Jdr x Jdr 2 {k d (r)[l + (P + 2)U( ri ,r 2 )/P} 

+n(r 1; r 2 ) + 6(n, r 2 ) + c(r h r 2 )}V K {r 1 )V K (r 2 ) ■ (15) 



The leading order contribution to (15) comes from il(r 1 ,r 2 ). Using the expansion of 



n(r 1 ,r 2 ) in terms of the diffusion modes 0q, the respective contribution is 

, 4 1 U dr4> Q (r)VJr)] 2 A 2 A 2 

a M = wv v — — = a i3f ■ ( 16) 

r 1 2 

For a disk in 2D the coefficient a = (2n)^ 1 Y.q I dr<f> q (r)(l — r 2 ) -1 ^ 2 /q 2 , where <p q 
are the diffusion modes in a disk of radius R — 1 (Q — R~ l q). The short-range 
term kd in (|T5|) contributes a term ~ (£ / L) 2 (A 2 / (5 g) that for a diffusive dot is para- 



metrically small compared with (jig) , but in the ballistic limit its contribution be- 



comes comparable to flip]). The surface charge contribution to the variance ([16]) dom- 
inates the bulk contribution ([8]). The bulk-surface cross- correlation term is given by 
(2AA/V 2 ) / dn J dr 2 \V 2 (\i) a (r 1 )\ 2 \i) a (r 2 )\ 2 ) - 1] V K (r 2 ). The dominating contribution due 
to kd (~ (£/L) 2 A 2 /g) is suppressed in a diffusive dot but is comparable to ([RJ) in the ballistic 
limit. 

Similarly, the covariance of v ai and v a g (due to surface charge) is evaluated to be 



^sjdnj dr 2 [(V 2 (\Mri)\ 2 \Mr2)\ 2 ) - l) + 3 (V 2 {^{r^^ 2 ) - l)] V K ( ri )V K (r 2 



(3V 2 



dr\ J dr 2 



k d (r) (1 + M±ln(n, r 2 ) J + n(n, r 2 ) + 6(n, r 2 ) + 3c(n, r 2 ) 



17) 



xK(n)K(r 2 ) . 



The leading order contributions for both a diffusive and ballistic dots are similar to the ones 
contributing to the variance, and comparing ([T7|) with (pT5|), we find 

- V^i « ^CT 2 (^a 7 ) • (18) 

We therefore obtain the following relation 

a 2 (5e^) = I m ( m + l)a 2 (5e«) = im(m + l)a^ . (19) 

Thus, for m ^> 1, the FSD of a level is linear in the number of added electrons m, 
in agreement with the result (|l|) of the parametric approach. The scrambling parameter 
8x\ can be determined by comparing the variance ([!]) in the parametric approach with the 
corresponding variance (|H|) in the microscopic approach 

M4f- (20) 

A complete scrambling of the spectrum corresponds to a parametric change of one avoided 
crossing, i.e., 5x m ~ 1. Using (pOj), we find 

m c ~ ~ (/^) 1/2 . (21) 



When an interaction scale A is introduced (e.g., screening due to external gates), we find 



In a recent experiment |?|] the standard deviation a p (m) = a(Vj^ +rn — VV), where Vj^ is 
the position in gate voltage for the jV-th peak, was measured as a function of the separation 
m in peak number. In Koopmans' limit V_\f +m — V^- ~ ^f+m 1 ~ e Af" > = ^m +m ^ + ^ e j^\ 
where A^"" 1 "" 1 ) = e^f+m 1 ~ e^ +m ' ) is the distance between levels separated by m consecutive 
spacings in a dot with a fixed number of electrons (Af+m). The latter quantity is unrelated to 
scrambling and its fluctuation for m < g is completely determined from RMT: cr(A,^ +m * ) ) oc 
A lnm. In principle one can extract a(5effi) from the measured a p (m). In practice this is 
difficult to do with the available experimental statistics since for m < m c , a(5effi) is smaller 
than cr(A^ +m )), and a p (m) is dominated by the lnm term of RMT. Nevertheless, a p (m), 
when plotted versus lnm, appears to follow a straight-line behavior with a small concave 
deviation ||25|| , in qualitative agreement with an additional contribution due to scrambling. 

In conclusion, we have studied spectral scrambling when several electrons are added to a 
Coulomb-blockade quantum dot. In the absence of surface charge the fluctuation standard 
deviation of a level upon the addition of m electrons behaves like y/m, but when surface 
charge is present this standard deviation is linear in m. The latter result is the one obtained 
in the conventional parametric random matrix approach. Our analysis was carried out 
within the framework of the HF-Koopmans picture and was limited to energy fluctuations. 
A microscopic description of wave functions scrambling is an outstanding problem. 
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